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Abstract 

Given two points in a simple polygon P of n vertices, its geodesic 
distance is the length of the shortest path that connects them among all 
paths that stay within P. The geodesic center of P is the unique point in 
P that minimizes the largest geodesic distance to all other points of P. In 
1989, Pollack, Sharir and Rote [Disc. & Comput. Geom. 89] showed an 
0(n log n)-time algorithm that computes the geodesic center of P. Since 
then, a longstanding question has been whether this running time can 
be improved (explicitly posed by Mitchell [Handbook of Computational 
Geometry, 2000]). In this paper we affirmatively answer this question and 
present a linear time algorithm to solve this problem. 


1 Introduction 

Let P be a simple polygon with n vertices. Given two points x,y in P, the 
geodesic path n(x, y ) is the shortest-path contained in P connecting x with y. If 
the straight-line segment connecting x with y is contained in P, then n(x,y) is 
a straight-line segment. Otherwise, tt( x,y) is a polygonal chain whose vertices 
(other than its endpoints) are reflex vertices of P. We refer the reader to !I9| 
for more information on geodesic paths refer. 

The geodesic distance between x and y, denoted by \-K(x,y)\, is the sum 
of the Euclidean lengths of each segment in Tr(x,y). Throughout this paper, 
when referring to the distance between two points in P, we refer to the geodesic 
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distance between them. Given a point x £ P, a (geodesic) farthest neighbor 
of x, is a point f P (x) (or simply f(x)) of P whose geodesic distance to x is 
maximized. To ease the description, we assume that each vertex of P has a 
unique farthest neighbor. We can make this general position assumption using 
simulation of simplicity [9]. 

Let F P (x) be the function that, for each x £ P, maps to the distance to 
a farthest neighbor of x (i.e., F P {x) = \n(x, f(x))\). A point cp £ P that 
minimizes F P {x) is called the geodesic center of P. Similarly, a point s £ P 
that maximizes F P {x) (together with its farthest neighbor) is called a geodesic 
diametral pair and their distance is known as the geodesic diameter. Asano and 
Toussaint [5] showed that the geodesic center is unique (whereas it is easy to 
see that several geodesic diametral pairs may exist). 

In this paper, we show how to compute the geodesic center of P in 0{n) 
time. 

1.1 Previous Work 

Since the early 1980s the problem of computing the geodesic center (and its 
counterpart, the geodesic diameter) has received a lot of attention from the 
computational geometry community. Chazelle [7] gave the first algorithm for 
computing the geodesic diameter (which runs in 0(n 2 ) time using linear space). 
Afterwards, Suri [24] reduced it to 0(nlogn)-time without increasing the space 
constraints. Finally, Hershberger and Suri m presented a fast matrix search 
technique, one application of which is a linear-time algorithm for computing the 
diameter. 

The first algorithm for computing the geodesic center was given by Asano 
and Toussaint [3], and runs in 0(n 4 log n)-time. In 1989, Pollack, Sharir, and 
Rote [22] improved it to 0(n log n) time. Since then, it has been an open 
problem whether the geodesic center can be computed in linear time (indeed, 
this problem was explicitly posed by Mitchell [f§i Chapter 27]). 

Several other variations of these two problems have been considered. Indeed, 
the same problem has been studied under different metrics. Namely, the L\ 
geodesic distance [6j , the link distance (231 El II] (where we look for the path 
with the minimum possible number of bends or links), or even rectilinear link 
distance [20J[2T| (a variation of the link distance in which only isothetic segments 
are allowed). The diameter and center of a simple polygon for both the L\ and 
rectilinear link metrics can be computed in linear time (whereas 0(n log n) time 
is needed for the link distance). 

Another natural extension is the computation of the diameter and center 
in polygonal domains (i.e., polygons with one or more holes). Polynomial time 
algorithms are known for both the diameter [5] and center [5|, although the 
running times are significantly larger (i.e., 0(n 7 73 ) and 0(n 12+s ), respectively). 
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1.2 Outline 


In order to compute the geodesic center, Pollack et al. [22j introduce a linear 
time chord-oracle. Given a chord C that splits P into two sub-polygons, the 
oracle determines which sub-polygon contains cp. Combining this operation 
with an efficient search on a triangulation of P, Pollack et al. narrow the search 
of cp within a triangle (and find the center using optimization techniques). 
Their approach however, does not allow them to reduce the complexity of the 
problem in each iteration, and hence it runs in 0(nlogn) time. 

The general approach of our algorithm described in Section [6] is similar: 
partition P into 0(1) cells, use an oracle to determine which cell contains cp, 
and recurse within the cell. Our approach differs however in two important 
aspects that allows us to speed-up the algorithm. First, we do not use the 
chords of a triangulation of P to partition the problem into cells. We use 
instead a cutting of a suitable set of chords. Secondly, we compute a set $ of 
0{n) functions, each defined in a triangular domain contained in P, such that 
their upper envelope, <fi(x), coincides with F P {x). Thus, we can “ignore” the 
polygon P and focus only on finding the minimum of the function <p(x). 

The search itself uses e-nets and cutting techniques, which certify that both 
the size of the cell containing cp and the number of functions of <f> defined 
in it decrease by a constant fraction (and thus leads to an overall linear time 
algorithm). This search has however two stopping conditions, (1) reach a sub¬ 
problem of constant size, or (2) find a triangle containing cp. In the latter 
case, we show that <f>(x) is a convex function when restricted to this triangle. 
Thus, finding its minimum becomes an optimization problem that we solve in 
Section [ 7 ] using cuttings in R 3 . 

The key of this approach lies in the computation of the functions of $ and 
their triangular domains. Each function g(x) of $ is defined in a triangular 
domain A contained in P and is associated to a particular vertex w of P. 
Intuitively speaking, g{x ) maps points in A to their (geodesic) distance to w. 
We guarantee that, for each point x £ P, there is one function g defined in a 
triangle containing x , such that g{x) = F P {x). To compute these triangles and 
their corresponding functions, we proceed as follows. 

In Section [3j we use the matrix search technique introduced by Hershberger 
and Suri tl3) to decompose the boundary of P, denoted by dP, into connected 
edge disjoint chains. Each chain is defined by either (1) a consecutive list of 
vertices that have the same farthest neighbor v (we say that v is marked if it 
has such a chain associated to it), or (2) an edge whose endpoints have different 
farthest neighbors (such edge is called a transition edge). 

In Section |4j we consider each transition edge ab of dP independently and 
compute its hourglass. Intuitively, the hourglass of ab, H a b, is the region of P 
between two chains, the edge ab and the chain of dP that contains the farthest 
neighbors of all points in ab. Inspired by a result of Suri [2'4j . we show that 
the sum of the complexities of each hourglass defined on a transition edge is 
0{n). In addition, we provide a new technique to compute all these hourglasses 
in linear time. 
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Figure 1: Given two edge disjoint transition chains, their hourglasses are open and 
the bottom chains of their hourglasses are also edge disjoint. Moreover, these bottom 
chains appear in the same cyclic order as the top chains along dP. 


In Section [5] we show how to compute the functions in $ and their respective 
triangles. We distinguish two cases: (1) Inside each hourglass H a b of a transition 
edge, we use a technique introduced by Aronov et al. [2j that uses the shortest- 
path trees of a and b in H a b to decompose H a b into 0(\H a b\) triangles with their 
respective functions (for more information on shortest-path trees refer to [ID]). 
(2) For each marked vertex v we compute triangles that encode the distance 
from v. Moreover, we guarantee that these triangles cover every point of P 
whose farthest neighbor is v. Overall, we compute the 0{n) functions of $ in 
linear time. 


2 Hourglasses and Funnels 

In this section, we introduce the main tools that are going to be used by the 
algorithm. Some of the results presented in this section have been shown before 
in different papers. 

2.1 Hourglasses 

Given two points x and y on dP, let dP{x, y ) be the polygonal chain that starts 
at x and follows the boundary of P clockwise until reaching y. 

For any polygonal chain C = dP(po,pi,... ,p k ), the hourglass of C, de¬ 
noted by He, is the simple polygon contained in P bounded by C, ir{pk, f(po)), 
dP(fipo), f(pk)) and tt {f{p k ),Po)\ see Figure []] We call C and dP{f(p 0 ),f(p k )) 
the top and bottom chains of He, respectively, while 7r (p k , f(po)) and 7 T(f(p k ),Po) 
are referred to as the walls of He- 

We say that the hourglass He is open if its walls are vertex disjoint. We say 
C is a transition chain if f(po) 7 ^ f(p k ) and neither f(po) nor f(pk) are interior 
vertices of C. In particular, if an edge ab of dP is a transition chain, we say 
that it is a transition edge (see Figure [I]). 
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Lemma 2.1. [Rephrase of Lemma 3.1.3 of fflf] If C is a transition chain of 
dP, then the hourglass He is an open hourglass. 

Note that by Lemma |2.1[ the hourglass of each transition chain is open. 
In the remainder of the paper, all the hourglasses considered are defined by a 
transition chain, i.e., they are open and their top and bottom chains are edge 
disjoint. 

The following lemma is depicted in Figure [l] and is a direct consequence of 
the Ordering Lemma proved by Aronov et al. [2j Corollary 2.7.4]. 

Lemma 2 . 2 . Let C\, O 2 , O 3 be three edge disjoint transition chains of dP that 
appear in this order when traversing clockwise the boundary of P. Then, the 
bottom chains of Hc 1 , Hc 2 and Hq 3 are also edge disjoint and appear in this 
order when traversing clockwise the boundary of P. 

Let 7 be a geodesic path joining two points on the boundary of P. We say 
that 7 separates two points 27 and x 2 of dP if the points of X = { 27 , X 2 } and 
the endpoints of 7 alternate along the boundary of P (27 and 27 could coincide 
with the endpoints of 7 in degenerate cases). We say that a geodesic path 7 
separates an hourglass H if it separates the points of its top chain from those 
of its bottom chain. 

Lemma 2.3. Let Cj,, C r be edge disjoint transition chains of dP. Then, 
there is a set of t = 0(1) geodesic paths 71 ,. .., 7 t with endpoints on dP such 
that for each 1 < i < r there exists 1 < j < t such that 7 j separates Hci ■ 
Moreover, this set can be computed in 0(n) time. 

Proof. Aronov et al. showed that there exist four vertices 17,...,17 of P and 
geodesic paths 27 ), n(v 2 , 03 ), tt(v 3 , 17 ) such that for any point x £ dP, one 
of these paths separates x from f{x) [U Lemma 2.7.6]. Moreover, they show 
how to compute this set in O(ro) time. 

Let F = {n(vi, Vj) : 1 < * < j < 4} and note that V\,..., 17 split the bound¬ 
ary of P into at most four connected components. If a chain Ci is completely 
contained in one of these components, then one path of T separates the top and 
bottom chain of He, ■ Otherwise, some vertex Vj is an interior vertex of Ci. 
However, because the chains C\,..., C r are edge disjoint, there are at most four 
chains in this situation. For each chain Ci containing a vertex Vj, we add the 
geodesic path connecting the endpoints of Ci to T. Therefore, T consists of 0(1) 
geodesic paths and each hourglass He, has its top and bottom chain separated 
by some path of T. Since only 0(1) additional paths are computed, this can be 
done in linear time. □ 

A chord of P is an edge joining two non-adjacent vertices a and b of P such 
that ab C P. Therefore, a chord splits P into two sub-polygons. 

Lemma 2.4. [Rephrase of Lemma 3-4-3 of [ffj] Let C \,..., C r be a set of edge 
disjoint transition chains of dP that appear in this order when traversing clock¬ 
wise the boundary of P. Then each chord of P appears in 0(1) hourglasses 
among Hc 1 , • • •, Hc r ■ 
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Proof. Note that chords can only appear on walls of hourglasses. Because hour¬ 
glasses are open, any chord must be an edge on exactly one wall of each of these 
hourglasses. Assume, for the sake of contradiction, that there exists two points 
s,t £ P whose chord st is in three hourglasses Hc i ,Hc j and Hc k (for some 
1 <i<j<k<r) such that s visited before t when going from the top chain 
to the bottom one along the walls of the three hourglasses. Let s,; and U be 
the points in the in the top and bottom chains of Hci , respectively, such that 
7 r(si,ti) is the wall of He, that contains st (analogously, we define Sk and tk) 
Because Cj lies in between C( and Ck , Lemma 2.2 implies that the bottom 
chain of Cj appears between the bottom chains of C, and Ck ■ Therefore, Cj 
lies between Si and Sk and the bottom chain of Hcj lies between ti and tk- 
That is, for each x £ Cj and each y in the bottom chain of He ., the geodesic 
path n(x,y) is “sandwiched” by the paths 7r (si,U) and ir(sk,tk )■ In particular, 
n{x,y) contains st for each pair of points in the top and bottom chain of Hc r 
However, this implies that the hourglass He is not open—a contradiction that 
comes from assuming that st lies in the wall of three open hourglasses, when this 
wall is traversed from the top chain to the bottom chain. Analogous arguments 
can be used to bound the total number of walls that contain the edge st (when 
traversed in any direction) to 0(1). □ 


Lemma 2.5. Let x,u,y,v be four vertices of P that appear in this cyclic order 
in a clockwise traversal of dP. Given the shortest-path trees T x and T y of x and 
y in P, respectively, such that T x and T y can answer lowest common ancestor 
(LCA) queries in 0(1) time, we can compute the path tt { u , v ) in 0(|7r(zt, u)|) 
time. Moreover, all edges ofir{u,v), except perhaps one, belong to T x LiT y . 

Proof. Let X {resp. Y) be the set containing the LCA in T x ( resp. T y ) of u,y, 
and of v, y {resp. u, x and x, y). Note that the points of X U Y lie on the path 
Tt{x, y) and can be computed in 0(1) time by hypothesis. Moreover, using LCA 
queries, we can decide their order along the path n(x, y) when traversing it 
from x to y. (Both X and Y could consist of a single vertex in some degenerate 
situations). Two cases arise: 

Case 1. If there is a vertex x* £ X lying after a vertex y* £ F along tt { x , y), 
then the path 77(77, v) contains the path 7 r{y*, x*). In this case, the path 77(77, v ) 
is the concatenation of the paths 77(77,7/*), tt {y*,x*), and tt { x *, v ) and that the 
three paths are contained in T x U T y . Moreover, 77(77,77) can be computed in 
time proportional to its length by traversing along the corresponding tree; see 
Figure [2] (top). 

Case 2. In this case the vertices of X appear before the vertices of Y along 
77(37,7/). Let x' {resp. y') be the vertex of X {resp. Y) closest to x {resp. y). 

Let v! be the last vertex of ir{u,x) that is also in 77(77,7/). Note that v! 
can be constructed by walking from v! towards x y until the path towards y 
diverges. Thus, v! can be computed in 0 ( 177 ( 11 , 7 /)!) time. Define v' analogously 
and compute it in 0(|7r(n,t/)|) time. 

Let P' be the polygon bounded by the geodesic paths 77 ( 37 ', it'), 7 t { u ', 7 /'), 77 ( 7 /', v') 
and 77(77', x'). Because the vertices of X appear before those of Y along 77(3;, y), 

P’ is a simple polygon; see Figure [ 2 ] (bottom). 
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Figure 2: (top) Case 1 of the proof of Lemma 2.5 where the path n(u,v) contains 
a portion of the path n(x, y). (bottom) Case 2 of the proof of Lemma 2.5 where the 
path 7 r(u, n) has exactly one edge being the tangent of the paths 7r(u', y') and ir(v', x'). 


In this case the path n(u,y) is the union of w(u, u'), v') and n(v',v). 

Because tt(u,u') and tt(v',v) can be computed in time proportional to their 
length, it suffices to compute tt(u',v') in 0(\tt(u', t/)|) time. 

Note that P' is a simple polygon with only four convex vertices x', u', y' and 
v', which are connected by chains of reflex vertices. Thus, the shortest path 
from x' to y' can have at most one diagonal edge connecting distinct reflex 
chains of P'. Since the rest of the points in lie on the boundary of P' 

and from the fact that each edge of P' is an edge of T x U T y , we conclude all 
edges of 7r(w, u), except perhaps one, belong to T x UT y . 

We want to find the common tangent between the reflex paths ir(u ', x') and 
7t(t/, y'), or the common tangent of 7 t(u\ y') and 7t(t/, x’) as one of them belongs 
to the shortest path tt(u', v’). Assume that the desired tangent lies between the 
paths 7 t(m', x’) and n(v\ y'). Since these paths consist only of reflex vertices, the 
problem can be reduced to finding the common tangent of two convex polygons. 
By slightly modifying the linear time algorithm to compute this tangents, we 
can make it run in 0{\t:(u' ,v')\) time. 

Since we do not know if the tangent lies between the paths Tr(u',x') and 
7 T(v',y'), we process the chains Tr(u',y') and n(v',x') in parallel and stop when 
finding the desired tangent. Consequently, we can compute the path tt(u,v) in 
time proportional to its length. □ 

Lemma 2.6. Let P be a simple polygon with n vertices. Given k disjoint 
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transition chains C \,..., Ck of dP, it holds that 


El HcA = 0{n). 


Proof. Because the given transition chains are disjoint, Lemma 1 2.2 1 implies that 
the bottom chains of their respective hourglasses are also disjoint. Therefore, 
the sum of the complexities of all the top and bottom chains of these hourglasses 
is 0(n). To bound the complexity of their walls we use Lemma 2.4 Since no 


chord is used more than a constant number of times, it suffices to show that the 
total number of chords used by all these hourglasses is 0(n). 


To prove this, we use Lemma 2.3 to construct 0(1) split chains 71 ,... , 7 1 
such that for each 1 < * < fc, there is a split chain 7 j that separates the top and 
bottom chains of Hq • For each 1 < j <t, let 


W = {Hci : the top and bottom chain of He., are separated by 7 j}. 

Since the complexity of the shortest-path trees of the endpoints of 7 j is O(n) [ID], 
and from the fact that the chains Cj,.. •, G\ are disjoint, Lemma |2.5| implies 
that the total number of edges in all the hourglasses of W is 0(n). Moreover, 
because each of these edges appears in 0 ( 1 ) hourglasses among Ci,... ,Ck, we 
conclude that 

E \ H \ = °W- 

ffeW 

Since we have only 0(1) split chains, our result follows. O 


2.2 Funnels 

Let C = (po,... ,pk) be a chain of dP and let v be a vertex of P not in C. 
The funnel of v to C, denoted by S V (C), is the simple polygon bounded by 
O, n(pk,v) and n(v,po); see Figure [3] (a). Note that the paths n(v,pk) and 
n(v,p 0 ) may coincide for a while before splitting into disjoint chains. See Lee 
and Preparata [T5] or Guibas et al. m for more details on funnels. 

A subset R C P is geodesically convex if for every x,y £ R, the path n(x, y) 
is contained in R. This funnel S v (C) is also known as the geodesic convex hull 
of C and v, i.e., the minimum geodesically convex set that contains v and C. 

Given two points x,y £ P, the (geodesic) bisector of x and y is the set of 
points contained in P that are equidistant from x and y. This bisector is a curve, 
contained in P, that consists of circular arcs and hyperbolic arcs. Moreover, this 
curve intersects dP only at its endpoints [Tj Lemma 3.22]. 

The (farthest) Voronoi region of a vertex v of P is the set of points R(v) = 
{x £ P : F P (x) = | 7 r(x,u)|} (including boundary points). 

Lemma 2.7. Let v be a vertex of P and let C be a transition chain such 
R(v) D dP C C and v ^ C. Then, R(v) is contained in the funnel S V (C) 




Proof. Let a and b be the endpoints of C such that a, b, f(a) and f(b) appear 
in this order in a clockwise traversal of dP. Because R(v) D dP C C, we know 
that v lies between /(a) and f(b). 

Let a ( resp. (3) be the bisector of v and /(a) ( resp. f(b)). Let h a ( resp. hb) 
be the set of points of P that are farther from v than from /(a) (resp. f(b)). 
Note that a is the boundary of h a while (3 bounds hb- 

By definition, we know that R(v) C h a D hb- Therefore, it suffices to show 
that h a nh b C S V (C). Assume for a contradiction that there is a point of h a nhb 
lying outside of S„ (C). By continuity of the geodesic distance, the boundaries 
of h a fl hb and S V (C) must intersect. Because a ^ h a and b ^ h b , both bisectors 
a and (3 must have an endpoint on the edge ab. Since the boundaries of h a D h b 
and S„(C) intersect, we infer that (3 fl n(v, b) ^ 0 or a. fl n(v, a) ^ 0. Without 
loss of generality, assume that there is a point w £ (3 C I tt(v, b), the case where 
w lies in afl n(v, a) is analogous. 

Since w € (3, we know that | 7 r(w,u)| = \w(w, f(b))\. By the triangle inequal¬ 
ity and since w cannot be a vertex of P as w intersects dP only at its endpoints, 
we get that 

\n(b,f(b)) | < |tt(&,w)| + \n(w,f(b))\ = \n(b,w)\ + \n(w,v)\ = |tt(6,i;)|. 

Which implies that b is farther from v than from f(b) —a contradiction that 
comes from assuming that h a n h b is not contained in 5'„(C). □ 


3 Decomposing the boundary 

In this section, we decompose the boundary of P into consecutive vertices that 
share the same farthest neighbor and edges of P whose endpoints have distinct 
farthest neighbors. 

Using a result from Hershberger and Suri (13) . in 0(n) time we can compute 
the farthest neighbor of each vertex of P. Recall that the farthest neighbor of 
each vertex of P is always a convex vertex of P [3] and is unique by our general 
position assumption. 

We mark the vertices of P that are farthest neighbors of at least one vertex 
of P. Let M denote the set of marked vertices of P (clearly this set can be 
computed in 0(n ) time after applying the result of Hershberger and Suri). In 
other words, M contains all vertices of P whose Voronoi region contains at least 
one vertex of P. 

Given a vertex v of P, the vertices of P whose farthest neighbor is v appear 
contiguously along dP [2- Therefore, after computing all these farthest neigh¬ 
bors, we effectively split the boundary into subchains, each associated with a 
different vertex of M; see Figure [3] (b). 

Let a and b be the endpoints of a transition edge of dP such that a appears 
before b in the clockwise order along dP. Because ab is a transition edge, we 
know that /(a) ^ f(b). Recall that we have computed /(a) and f(b) in the 
previous step and note that f(a) appears also before f(b) along this clockwise 
order. For every vertex v that lies between f(a) and f(b) in the bottom chain 
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Figure 3: a) The funnel ^(C) of a vertex v and a chain C contained in dP are de¬ 
picted. b ) Each vertex of the boundary of P is assigned with a farthest neighbor which 
is then marked. The boundary is then decomposed into vertex disjoint chains, each 
associated with a marked vertex, joined by transition edges (blue) whose endpoints 
have different farthest neighbors. 


of i? a b, we know that there cannot be a vertex u of P such that f{u) = v. As 
proved by Aronov et al. [2] Corollary 2.7.4], if there is a point x on dP whose 
farthest neighbor is v, then x must lie on the open segment (a, b). In other 
words, the Voronoi region R(v) restricted to dP is contained in (a, b). 


4 Building hourglasses 


Let E be the set of transition edges of dP. Given a transition edge ab £ E, 
we say that H a i, is a transition hourglass. In order to construct the triangle 
cover of P, we construct the transition hourglass of each transition edge of E. 
By Lemma 2.6 we know that J2 a b£E \H a b\ = 0(n). Therefore, our aim is to 


compute the cover in time proportional to the size of H^. 


By Lemma 2.3 we can compute a set of 0(1) separating paths such that 


for each transition edge ab 7 the transition hourglass H a \, is separated by one 
(or more) paths in this set. For each endpoint of the 0(1) separating paths 
we compute its sliortest-path tree [lOj. In addition, we preprocess these trees 
in linear time to support LCA queries [12] . Both computations need linear 
time per endpoint and use 0(n) space. Since we do this process for a constant 
number of endpoints, overall this preprocessing takes O(n) time. 

Let 7 be a separating path whose endpoints are x and y. Note that 7 
separates the boundary of P into two chains S and S' such that SUS' = dP. Let 
T~L{ 7 ) be the set of each transition hourglass separated by 7 whose transition edge 
is contained in S (whenever an hourglass is separated by more than one path, 
we pick one arbitrarily). Note that we can classify all transition hourglasses into 
the sets H(y) in 0(n) time (since 0 ( 1 ) separating paths are considered). 

We claim that we can compute all transition hourglass of Hipf) in 0(n) time. 
By construction, the wall of each of these hourglasses consists of a (geodesic) 
path that connects a point in S with a point in S'. Let u £ S and v £ S' be 


10 








two vertices such that n(u,v) is the wall of a hourglass in R("f). Because LCA 
queries can be answered in 0(1) time [T2], Lemma 2.5 allows us to compute this 
path in 0(\tt(u,v)\) time. Therefore, we can compute all hourglasses of H{ 7) in 
\ H \ + n ) = 0(n) time by Lemma 


2.6 


Because only 0(1) separating 


paths are considered, we obtain the following result. 


Lemma 4.1. We can construct the transition hourglass of all transition edges 
of P in 0(n ) time. 


5 Covering the polygon with apexed triangles 

An apexed triangle A = (a, 6, c) with apex a is a triangle contained in P with 
an associated distance function g/^ix), called the apex function of A, such that 
(1) a is a vertex of P, (2) b, c £ dP , and (3) there is a vertex w of P, called the 
definer of A, such that 

, , _ J —00 if x A 

'^ A ( \xa\ + |7r(a, w;)| = \ir(x, u>)| if x £ A 

In this section, we show how to find a set of O(n) apexed triangles of P 
such that the upper envelope of their apex functions coincides with F P (x). To 
this end, we first decompose the transition hourglasses into apexed triangles 
that encode all the geodesic distance information inside them. For each marked 
vertex v € M we construct a funnel that contains the Voronoi region of v. We 
then decompose this funnel into apexed triangles that encode the distance from 
v. 

5.1 Inside the transition hourglass 

Let ab be a transition edge of P such that b is the clockwise neighbor of a along 
dP. Let B a b denote the bottom chain of H ab after removing its endpoints. As 
noticed above, a point on dP can be farthest from a vertex in B ab only if it lies 
in the open segment ab. That is, if v is a vertex of B ab such that R(v) ^ 0, then 
R(v ) fl dP C ab. 

In fact, not only this Voronoi region is inside H ab when restricted to the 
boundary of P, but also R{v) C H ab . The next result follows trivially from 
Lemma 12.71 

Corollary 5.1. Let v be a vertex of B ab . If R(v) ^ 0, then R(v) C H ab . 

Our objective is to compute 0{\H ab \) apexed triangles that cover H ab , each 
with its distance function, such that the upper envelope of these apex functions 
coincides with F P ( x ) restricted to H ab where it “matters”. 

The same approach was already used by Pollack et al. in [22] Section 3]. 
Given a segment contained in the interior of P, they show how to compute 
a linear number of apexed triangles such that F P {x ) coincides with the upper 
envelope of the corresponding apex functions in the given segment. 
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Figure 4: (left) A vertex v visible from the segment ab lying on the bottom chain of 
H a i,, and the triangle A„ which contains the portion of ab visible from v. (right) The 
children mi and M 2 of v are visible from ab while W 3 is not. The triangle A„ is split 
into apexed triangles by the rays going from ui and M 2 to v. 


While the construction we follow is analogous, we use it in the transition 
hourglass H ab instead of the full polygon P. Therefore, we have to specify what 
is the relation between the upper envelope of the computed functions and F P (x). 
We will show that the upper envelope of the apex functions computed in H ab 
coincides with F P (x) inside the Voronoi region R(v) of every vertex v £ B ab . 

Let T a and T b be the shortest-path trees in H ab from a and b , respectively. 
Assume that T a and T b are rooted at a and b, respectively. We can compute 
these trees in 0(\H ab \) time [TO] , For each vertex v between f(a) and /(&), let 
v a and v b be the neighbors of v in the paths n(v,a) and 7 r(v,b), respectively. 
We say that a vertex v is visible from ab if v a 7^ v b . Note that if a vertex is 
visible, then the extension of these segments must intersect the top segment 
ab. Therefore, for each visible vertex v, we obtain a triangle A„ as shown in 
Figure [fj 

We further split A t , into a series of triangles with apex at v as follows: Let 
it be a child of v in either T a or T b . As noted by Pollack et ah, v can be of three 
types, either (1) u is not visible from ab (and is hence a child of v in both T a 
and Tb); or (2) u is visible from ab, is a child of v only in T b , and v b vu is a left 
turn; or (3) u is visible from ab, is a child of v only in T a , and v a vu is a right 
turn. 

Let u\,... ,Uk- 1 be the children of v of type (2) sorted in clockwise order 
around v. Let c(v ) be the maximum distance from v to any invisible vertex in 
the subtrees of T a and T b rooted at v; if no such vertex exists, then c(v) = 0. 
Define a function d; ( v ) on each vertex v of H ab in a recursive fashion as follows: 
If v is invisible from ab, then di(v) = c(v). Otherwise, let di(v) be the maximum 
of c(v) and max{di(ui) + \uiv\ : Ui is a child of v of type (2)}. Similarly we define 
a symmetric function d r (v) using the children of type (3) of v. 

For each 1 < i < k — 1, extend the segment UiV past v until it intersects ab 
at a point Sj. Let so and be the intersections of the extensions of vv a and vv b 
with the segment ab. We define then k triangles contained in A„ as follows. For 
each 0 < i < k — 1, consider the triangle A(s,, v, Sj+i) whose associated apexed 
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(left) function is 


f \xv\+maxj >i {c(v),\vuj\+di(uj)} if x £ A(s*, v, s i+ i) 

\ —oo otherwise 


In a symmetric manner, we define a set of apexed triangles induced by the type 
(3) children of v and their respective apexed (right) functions. 

Let gi,... ,g r and Ai, ..., A r respectively be an enumeration of all the gen¬ 
erated apex functions and triangles such that g t is defined in the triangle Aj. 
Because each function is determined uniquely by a pair of adjacent vertices 
in T a or in Tb, and since these trees have 0(\H a b\) vertices, we conclude that 
r = 0(\H ab \). 

Note that for each 1 < i < r, the triangle A i has two vertices on the 
segment ab and a third vertex, say cq, called its apex such that for each x £ Ai, 
gi(x) = \n{x,vji)\ for some vertex Wi of H a b- We refer to Wi as the definer of 
A j. Intuitively, A, defines a portion of the geodesic distance function from Wi 
in a constant complexity region. 

Lemma 5.2. Given a transition edge ab of P , we can compute a set A a b of 
0(\H a b\) apexed triangles in 0(\H a b\) time with the property that for any point 
p £ P such that f(p) £ B a b, there is an apexed triangle A £ A a b with apex 
function g and definer equal to f(p) such that 

1. p £ A and 


2. g(p) = F P {p). 


Proof. Because p £ R(f(p)), Lemma 5.1 implies that p £ H a b . Consider the 
path 7r(p, f(jp)) and let v be the neighbor of p along this path. By construction 
of A a b, there is a triangle A £ A a b apexed at v with definer w that contains 
p. The apex function g(x) of A encodes the geodesic distance from x to w. 
Because F P (x) is the upper envelope of all the geodesic functions, we know that 
g{p) < F P (p). 

To prove the other inequality, note that if v = f(p), then trivially g(p) = 
\pv\ + |7r(u,u;)| > \pv\ = 17r {jp, f (p))| = F P {p). Otherwise, let z be the next 
vertex after v in the path ir(p, f(p)). Three cases arise: 

(a) If z is invisible from ab, then so is f(p) and hence, 


\n(p,f(p))\ = | pv\ + \tt(v,/{p))\ < | pv\ + c(v) < g(p). 


(6) If z is a child of type (2), then z plays the role of some child Uj of v in 
the notation used during the construction. In this case: 


Ha/(p))| = IH + \vz\ + \n(z,f{p))\ < | pv\ + \vuj\ +di(uj) < g(p). 

(c) If z is a child of type (3), then analogous arguments hold using the (right) 
distance d r . 

Therefore, regardless of the case F P {p) = \n(p,f(p))\ < g(p). 
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To bound the running time, note that the recursive functions di , d r and c 
can be computed in 0(\T a \ + |Tj|) time. Then, for each vertex visible from ab , 
we can process it in time proportional to its degree in T a and T b . Because the 
sum of the degrees of all vertices in T a and T b is 0(|T o | + \T b \) and from the fact 
that both \T a \ and \T b \ are 0(\H ab \), we conclude that the total running time 
to construct A a b is 0(\H ab \). □ 

In other words, Lemma |5.2| says that no information on farthest neighbors 
is lost if we only consider the functions in A a b within H ab .In the next section 
we use a similar approach to construct a set of apexed triangles (and their 
corresponding apex functions), so as to encode the distance from the vertices of 

M. 


5.2 Inside the funnels of marked vertices 


Recall that for each marked vertex v £ M, we know at least of one vertex on dP 
such that v is its farthest neighbor. For any marked vertex v, let Ui,... ,Uk- 1 
be the vertices of P such that v = f(uf) and assume that they appear in this 
order when traversing dP clockwise. Let uq and Uk be the neighbors of u± and 
Uk -1 other than 112 and Uk- 2 , respectively. Note that both uqUi and Uk-iUk 
are transition edges of P. Thus, we can assume that their transition hourglasses 
have been computed. 

Let C v = (u 0 ,... ,Uk) and consider the funnel S V (C V ). We call C v the 
main chain of S V (C V ) while Tr(uk,v) and ir(v,uo) are referred to as the walls 
of the funnel. Because v = f(u\) = f(uk- 1 ), we know that v is a vertex 
of both H UoUl and H Uk lUk . By definition, we have tt(v,Uq) C H UqUi and 
7 r(v,Uk) C H Uk lUk . Thus, we can explicitly compute both paths 7r(v,uo) and 
tf (v,u k ) in O(\H U0Ul I + | H Uk 
constructed in 0(k + \H Uf)Ul \ + \H Uk _ lUk \) time. 

the total sum of the complexities of the transition hourglasses is 0(n 


. |) time. So, overall, the funnel S' t ,(C'„) can be 


Recall that, by Lemma 2.6 

hi 


particular, we can bound the total time needed to construct the funnels of all 
marked vertices by 0{n). 

Since the complexity of the walls of these funnels is bounded by the com¬ 
plexity of the transition hourglasses used to compute them, we get that 

\S v (C v )\ = o(n+ J2 =0(n). 

v&M \ ab£E ) 

Lemma 5.3. Let x be a point in P. If v = f{x) is a vertex of M, then 
x £ S„{C V ). 

Proof. Since /(mq) 7 ^ f{uk), C v is a transition chain. Moreover, C v contains 


R(v) H dP by definition. Therefore, Lemma 2.7 implies that R(v) C S V (C V ). 
Since v = f{x), we know that x £ R(v) and hence that x £ S„(C V ). □ 

We now proceed to split a given funnel into 0(|5„(C„)|) apexed triangles that 
encode the distance function from v. To this end, we compute the shortest-path 
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tree T v of v in S„(C V ) in 0(|S'„(C'„)|) time [IT]. We consider the tree T v to be 
rooted at v and assume that for each node u of this tree we have stored the 
geodesic distance |7r(zt,u)|. 

Start an Eulerian tour from v walking in a clockwise order of the edges. Let 
Let w\ be the first leaf of T v found, and let W 2 and W 3 be the next two vertices 
visited in the traversal. Two cases arise: 

Case 1: w±,W 2 ,W 3 makes a right turn. We define s as the first point hit by 
the ray apexed at w2 that shoots in the direction opposite to W3. 

We claim that w 1 and s lie on the same edge of the boundary of S^(C V ). 
Otherwise, there would be a vertex u visible from w2 inside the wedge with 
apex w2 spanned by w 1 and W3. Note that the first edge of the path n(u,v) 
is the edge UW2- Therefore, uw2 belongs to the shortest-path T v contradicting 
the Eulerian order in which the vertices of this tree are visited as u should be 
visited before W3. Thus, s and w\ lie on the same edge and s can be computed 
in 0 ( 1 ) time. 

At this point, we construct the apexed triangle A{w2,w\, s) apexed at W2 
with apex function 

, f \xu>2\ + \'k(w2,v)\ if X G A(W2,Wi,s) 

' ( —00 otherwise 

We modify tree T v by removing the edge W 1 W 2 and replacing the edge W 3 W 2 by 
the edge W3S; see Figure [5] 

Case 2: W\,W 2 ,W 3 makes a left turn and w\ and W 3 are adjacent, then 
if w 1 and W3 lie on the same edge of dP , we construct an apexed triangle 
A(w2,w± 1 W3) apexed at w2 with apex function 

(x)= ( \xw 2 \ + \tt(w 2 ,v)\ if X G A(w2,Wx,W3) 

\ —00 otherwise 


Otherwise, let s be the first point of the boundary of S V (C V ) hit by the ray 
shooting from W 3 in the direction opposite to w 2 . 

By the same argument as above, we can show that w\ and s lie on the same 
edge of the boundary of S„(C V ) (and thus, we can compute s in 0(1) time). We 
construct an apexed triangle A(u> 2 , w'i, s) apexed at w 2 with apex function 

/ x_ f \xu> 2 \ + \'k(w 2 ,v)\ if x € A(w 2 ,W!, s) 
y —00 otherwise 

We modify the tree T v by removing the edge W 1 W 2 and adding the edge W 3 S; 
see Figure [5] for an illustration. 

Lemma 5.4. The above procedure runs in 0(|S'„(C'„)|) time and computes 
0(|5'„(C'^)|) interior disjoint apexed triangles such that their union covers S^(C V ). 
Moreover, for each point x € R(v), there is an apexed triangle A with apex func¬ 
tion g(x) such that ( 1 ) x G A and (2) g{x) = F P (x). 
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Figure 5: The funnel S V (C V ) and the shortest-path tree from v are depicted in (a) 
The two cases of the algorithm described in Lemma 5.4 are shown in (6) and (c). 


Proof. The above procedure splits S V (C V ) into apexed triangles, such that their 
apex function in each of them is defined as the geodesic distance to v. By 
Lemma 5.3 if x £ R{v), then x £ S„(C V ). Therefore, there is an apexed triangle 
A with apex function g(x) such that x £ A and g(x) = \n(x,v)\ = F P (x). 
Consequently, we obtain properties (1) and (2). 

We now bound the running time of the algorithm. The shortest-path tree 
T v from v is computed in 0(|S'„(C'„)|) time [TU]. For each leaf of T v we need 
a constant number of operations to determine in which of the cases we are in 
(and to treat it as well). Therefore, it suffices to bound the number of times 
these steps are performed. Note that a leaf is removed from the tree in each 
iteration. Since the number of leaves strictly decreases each time we are in Case 
2, this step cannot happen more than 0(|S'„(C u )|) times. In Case 1 a new leaf 
is added if Wi and W 3 do not lie on the same edge of dP. However, the number 
of leaves that can be added throughout is at most the number of edges of T v . 
Note that the edges added by either Case 1 or 2 are chords of the polygon and 
hence do not generate further leaves. Because \T V \ = 0(|S’„(C'„)|), we conclude 
that both Case 1 and 2 are only executed 0(|S„(C'„)|) times. □ 


6 Prune and search 

With the tools introduced in the previous sections, we can proceed to give 
the prune and search algorithm to compute the geodesic center. The idea of 
the algorithm is to partition P into 0(1) cells, determine on which cell of P the 
center lies and recurse on that cell as a new subproblem with smaller complexity. 

Naturally, we can discard all apexed triangles that do not intersect the new 
cell containing the center. Using the properties of the cutting, we can show that 
both the complexity of the cell containing the center, and the number of apexed 
triangles that intersect it decrease by a constant fraction in each iteration of 
the algorithm. This process is then repeated until either of the two objects has 
constant descriptive size. 
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Let r be the set all apexed triangles computed in previous sections. Lem¬ 


mas 


2.6 and 5.4 directly provide a bound on the complexity of r. 


Corollary 6.1. The set r consists of 0{n) apexed triangles. 

Let <j)(x ) be the upper envelope of the apex functions of every triangle in r 
(i.e., 4>{x) = m&x{gi(x) : gt(x) £ r}). The following result is a direct conse¬ 
quence of Lemmas 5.2 and |5.4[ and shows that the 0(n) apexed triangles of r 
not only cover P , but their apex functions suffice to reconstruct the function 
F P ( x). 

Lemma 6.2. The functions (j>(x) and F P {x) coincide in the domain of points 
of P, i.e., for each p £ P, (j>{p) = F P (p). 


Given a chord CofPa half-polygon of P is one of the two simple polygons in 
which C splits P. A f-cell of P is a simple polygon obtained as the intersection 
of at most four half-polygons. Because a 4-cell is the intersection of geodesically 
convex sets, it is also geodesically convex. 

Let R be a 4-cell of P and let tr be the set of apexed triangles of r that 
intersect R. Let m R = max{|i?|, |tr|}. Recall that, by construction of the 
apexed triangles, for each triangle of tr at least one and at most two of its 
boundary segments is a chord of P Let C be the set containing all chords that 
belong to the boundary of a triangle of tr. Therefore, |tr| < \C\ < 2|tr |. 

To construct an e-net of C, we need some definitions (for more information 
on e-nets refer to HI)- Let ip be the set of all open 4-cells of P. For each t £ ip, 
let C t = {C £ C : C fl t ^ 0} be the set of chords of C induced by t. Finally, let 
ipc = {Ct : t £ ip} be the family of subsets of C induced by <p. 

Let e > 0 (the exact value of e will be specified later). Consider the range 
space ( C,ipc ) defined by C and ipc ■ Because the VC-dimension of this range space 
is finite, we can compute an e-net N of (C, ipc ) in 0{n/e) = 0{n ) time jT7]. The 
size of AT is 0( j log {) = 0(1) and its main property is that any 4-cell that does 
not intersect a chord of N will intersect at most e\C\ chords of C. 

Observe that N partitions R into 0(1) sub-polygons (not necessarily 4- 
cells). We further refine this partition by performing a 4-cell decomposition. 
That is, we shoot vertical rays up and down from each endpoint of N, and 
from the intersection point of any two segments of N, see Figure [6] Overall, 
this partitions R into 0(1) 4-cells such that each either (?) is a convex polygon 
contained in P of at most four vertices, or otherwise (??) contains some chain 
of dP. Since \N\ = 0(1), the whole decomposition can be computed in 0{m R ) 
time (the intersections between segments of N are done in constant time, and 
for the ray shooting operations we walk along the boundary of R once). 

In order to determine which 4-cell contains the geodesic center of P, we 
extend each edge of a 4-cell to a chord C. This can be done with two ray¬ 
shooting queries (each of which takes 0{m R ) time). We then use the chord- 
oracle from Pollack et al. [25J Section 3] to decide which side of C contains cp. 
The only requirement of this technique is that the function F P (x) coincides with 
the upper envelope of the apex functions when restricted to C. Which is true 
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Figure 6: The e-net TV splits P into 0(1) sub-polygons that are further refined into a 
4-cell decomposition using 0(1) ray-shooting queries from the vertices of the arrange¬ 
ment defined by TV. 


by Lemma |6.2| and from the fact that Tp consists of all the apexed triangles of 
r that intersect R. 

Because the chord-oracle described by Pollack et al. [321 Section 3] runs in 
linear time on the number of functions defined on C, we can decide in total 
0(m R ) time on which side of C the geodesic center of P lies. Since our decom¬ 
position into 4-cells has constant complexity, we need to perform 0(1) calls to 
the oracle before determining the 4-cell R' that contains the geodesic center of 
P. 

The chord-oracle computes the minimum of F P (x) restricted to the chord 
before determining the side containing the minimum. In particular, if cp lies 
on any chord bounding R' , then the chord-oracle will find it. Therefore, we can 
assume that cp lies in the interior of R'. Moreover, since TV is a e-net, we know 
that at most e\C\ chords of C will intersect R'. 

Using a similar argument, we can show that the complexity of R' also de¬ 
creases: since \C\ < 2 \ tr \ < 2 m R , we guarantee that at most 2em R apexed 
triangles intersect R'. Moreover, each vertex of R' is in at least one apexed tri¬ 
angle of Tp by Lemma [b~2j and by construction, each apexed triangle can cover 
at most three vertices. Thus, by the pigeonhole principle we conclude that R' 
can have at most 6 em R vertices. Thus, if we choose e = 1/12, we guarantee 
that both the size of the 4-cell R' and the number of apexed triangles in Tp, are 
at most m R /2 . 

In order to proceed with the algorithm on R' recursively, we need to compute 
the set rp, with the at most e\C\ apexed triangles of Tp that intersect R' (i.e., 
prune the apexed triangles that do not intersect with R'). For each apexed 
triangle A £ Tp, we can determine in constant time if it intersects R' (either one 
of the endpoints is in R' C\dP or the two boundaries have non-empty intersection 
in the interior of P). Overall, we need 0(m R ) time to compute the at most e\C\ 
triangles of Tp that intersect R'. 

By recursing on R', we guarantee that after 0(logm R ) iterations, we reduce 
the size of either Tp or R' to constant. In the former case, the minimum of 
F P (x) can be found by explicitly constructing function <j> in 0(1) time. In the 
latter case, we triangulate R' and apply the chord-oracle to determine which 
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triangle will contain cp. The details needed to find the minimum of (j>(x) inside 
this triangle are giving the next section. 

Lemma 6.3. In 0{n ) time we can find either the geodesic center of P or a 
triangle containing the geodesic center. 


7 Solving the problem restricted to a triangle 

In order to complete the algorithm it remains to show how to find the geodesic 
center of P for the case in which R' is a triangle. If this triangle is in the interior 
of P, it may happen that several apexed triangles of r fully contain R' . Thus, 
the pruning technique used in the previous section cannot be further applied. 
We solve this case with a different approach. 

Recall that denotes the upper envelope of the apex functions of the 
triangles in r, and the geodesic center is the point that minimizes <p. The key 
observation is that, as it happened with chords, the function <j>(x) restricted to 
R' is convex. 

Let Ai,A 2 ,...,A m be the set of m = O(n) apexed triangles of r that 
intersect R'. Let g%{x) be the apex function of Aj such that 

, \ f \xaf\ + Ki if x £ A* 

9{x) = { -00 otherwise ’ 


where aj and Wi are the apex and the definer of Aj, respectively, and k, = 
\n(ai,Wi)\ is a constant. 

By Lemma 6.2 cf>(x) = F P (x). Therefore, the problem of finding the center 
is equivalent to the following optimization problem in R 3 : 

(PI). Find a point (x,r) £ R 3 minimizing r subject to x £ R' and 


gi(x) < r, for 1 < * < to. 

Thus, we need only to find the solution to (PI) to find the geodesic center 
of P. We use some remarks described by Megiddo in order to simplify the 
description of (PI) [181 . 

To simplify the formulas, we square the equation |xaj| < r — Ki'. 

||a ;|| 2 - 2x ■ at + ||aj || 2 = |xaj | 2 < (r - nf) 2 = r 2 - 2 r«j + k 2 . 

And finally for each 1 < * < to, we define the function /ij(x,r) as follows: 

7 / x / ||x|| 2 - 2x • di + ||a*|| 2 - r 2 + 2r«* - if x G A* 

—oo otherwise 


Therefore, our optimization problem can be reformulated as: 

(P2). Find a point (x, r) £ R 3 such that r is minimized subject to x £ R' 


hi(x , r) < 0 and r > max{Kj}, for 1 < i < m. 
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Let /i'(x,r) = \\x\\ 2 -2x-ai + \\ai\\ 2 —r 2 +2rm — nf be a function defined in the 
entire plane and let (P2') be an optimization problem analogous to (P2) where 
every instance of hi(x,r) is replaced by h[(x,r). The optimization (P2 7 ) was 
studied by Megiddo in [DJI . We provide some of the intuition used by Megiddo 
to solve this problem. 

Although the functions are not linear inside A,, they all have the 

same non-linear terms. Therefore, for i ^ j, we get that h[(x,r) = h'j(x,r) 
defines a separating plane 

7 i,j = {(ar,r) G R 3 : 2 {m - Kj)r - 2(a* - a,-) ■ x + ||a,i|| 2 - 11ay 11 2 - + k? = 0} 

As noted by Megiddo (18j , this separating plane has the following property: 
If the solution (x,r) to (P2') is known to lie to one side of 'ytj, then we know 
that one of the constraints is redundant. 

Thus, to solve (P2 7 ) it sufficed to have a side-decision oracle to determine on 
which side of a plane 7 ij the solution lies. Megiddo showed how to implement 
this oracle in a way that the running time is proportional to the number of 
constraints [15] , 

Once we have such an oracle, Megiddo’s problem can be solved using a prune 
and search similar to that introduced in Section [ 6 ] pair the functions arbitrarily, 
and consider the set of m/2 separating planes defined by these pairs. For some 
constant r, compute a 1/r-cutting in R 3 of the separating planes. A 1/r-cutting 
is a partition of the plane into 0 (r 2 ) convex regions each of which is of constant 
size and intersects at most m/2r separating planes. A cutting of planes can 
be computed in linear time in R 3 for any r = 0(1) Qj3]. After computing the 
cutting, determine in which of the regions the minimum lies by performing 0 ( 1 ) 
calls to the side-decision oracle. Because at least ( 1 — l)m/2r separating planes 
do not intersect this constant size region, for each of them we can discard one of 
the constraints as it becomes redundant. Repeating this algorithm recursively 
we obtain a linear running time. 

To solve (P2) we follow a similar approach, but our set of separating planes 
needs to be extended in order to handle apex functions as they are only defined 
in the same way as in (P2 7 ) in a triangular domain. Note that the vertices of 
each apexed triangle that intersect R! have their endpoints either outside of R' 
or on its boundary. 

7.1 Optimization problem in a convex domain 

In this section we describe our algorithm to solve the optimization problem 
(P2). To this end, we pair the apexed triangles arbitrarily to obtain m/2 pairs. 
By identifying the plane where P lies with the plane Z 0 = {(x,y,z) : z = 0}, 
we can embed each apexed triangle in R 3 . A plane-set is a set consisting of at 
most five planes in R 3 . For each pair of apexed triangles (A*, A j) we define 
a plane-set as follows: For each chord bounding either A, or A j, consider the 
line extending this chord and the vertical extrusion of this line in R 3 , i.e., the 
plane containing this chord orthogonal to Zq. Moreover, consider the separating 
plane 7 ij. The set containing these planes is the plane-set of the pair (A*, A j). 
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Let r be the union of all the plane-sets defined by the to/ 2 pairs of apexed 
triangles. Thus, T is a set that consists of 0(m) planes. Compute an 1/r-cutting 
of T in 0(m) time for some constant r to be specified later. Because r is constant, 
this 1 /r-cutting splits the space into 0 ( 1 ) convex regions, each bounded by a 
constant number of planes nn. Using a side-decision algorithm (to be specified 
later), we can determine the region Q of the cutting that contains the solution 
to (P2). Because Q is the region of a 1/r-cutting of T, we know that at most 
|T|/r planes of T intersect Q. In particular, at most |T|/r plane-sets intersect 
Q and hence, at least (r — l)|T|/r plane-sets do not intersect Q. 

Let (Aj, Aj) be a pair such that its plane-set does not intersect Q. Let Q' 
be the projection of Q on the plane Zq. Because the plane-set of this pair does 
not intersect Q, we know that Q' intersects neither the boundary of Aj nor that 
of Aj. Two cases arise: 

Case 1 . If either A j or Aj does not intersect Q' , then we know that their 
apex function is redundant and we can drop the constraint associated with this 
apexed triangle. 

Case 2. If Q' C A* f~l Aj, then we need to decide which constrain to drop. 
To this end, we consider the separating plane 7 jj. Notice that inside the vertical 
extrusion of Aj (~l Aj (and hence in Q ), the plane 7 jj has the property that if 
we know its side containing the solution, then one of the constraints can be 
dropped. Since 7 ^ does not intersect Q as 7 belongs to the plane-set of 
(Aj, Aj), we can decide which side of 7 contains the solution to (P2) and 
drop one of the constraints. 

Regardless of the case if the plane-set of a pair (Aj, A j) does not intersect 
Q, then we can drop one of its constraints. Since at least (r — l)|r|/r plane- 
sets do not intersect Q , we can drop at least (r — l)|r|/r constraints. Because 
|r| > m /2 as each plane-set contains at least one plane, by choosing r = 2 , we 
are able to drop at least |P|/2 > to/ 4 constraints. Consequently, after 0(m) 
time, we are able to drop to/ 4 apexed triangles. By repeating this process 
recursively, we end up with a constant size problem in which we can compute 
the upper envelope of the functions explicitly and find the solution to (P2) using 
exhaustive search. Thus, the running time of this algorithm is bounded by the 
recurrence T(m) = T(3m/4) + 0(m) which solves to 0(m). Because to = O(n), 
we can find the solution to (P2) in 0(n ) time. 

The last detail is the implementation of the side-decision algorithm. Given a 
plane 7 , we want to decide on which side lies the solution to (P2). To this end, 
we solve (P2) restricted to 7 , i.e., with the additional constraint of (x,r) £ 7 . 
This approach was used by Megiddo [IS], the idea is to recurse by reducing the 
dimension of the problem. Another approach is to use a slight modification of 
the chord-oracle described by Pollack et al. [22j Section 3]. 

Once the solution to (P2) restricted to 7 is known, we can follow the same 
idea used by Megiddo [T5| to find the side of 7 containing the global solution to 
(P2). Intuitively, we find the apex functions that define the minimum restricted 
to 7 . Since <j>{x) = F P (x) is locally defined by this functions, we can decide on 
which side the minimum lies using convexity. We obtain the following result. 


21 


Lemma 7.1 . Let R' be a convex trapezoid contained in P such that R' contains 
the geodesic center of P. Given the set of all apexed triangles of r that intersect 
R', we can compute the geodesic center of P in 0(n) time. 

The following theorem summarizes the result presented in this paper. 

Theorem 7.2. We can compute the geodesic center of any simple polygon P 
of n vertices in 0(n) time. 
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